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ABSTRACT 


The  non-linear  transformations  relating  rectangular 
coordinates  to  range  instrumentation  coordinates  for  con¬ 
ventional  radar,  phased  array  radar,  AME  (Angle  Measuring 
Equipment)  etc.,  are  considered.  The  transformations  are 
linearly  related  through  Jacobian  matrices  at  the  deri¬ 
vative  (rate)  and  first  differential  or  discrete  difference 
(residuals  or  linearizing  about  a  nominal)  level.  State 
variables  near  the  singularity  regions  of  these  matrices 
can  cause  numerical  computer  problems  for  the  unwary 
analyst. 
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INTRODUCTION 


The  development  of  operational  computer  programs  from 
math-ware  to  software  presents  many  opportunities  for  bad 
output  and  expensive  reruns.  Systematic  procedures  for 
obtaining  computable  math-ware  from  the  physical  problems 
of  range  testing  are  being  developed  in  the  state-space 
frame-work.  A  deeper  understanding  of  the  math-ware  will 
enhance  the  analysts'  ability  to  make  judgments  on  the 
practical  implications  of  the  quantitative  answers  pro¬ 
vided  by  the  digital  computer. 

Range  flight  testing  and  trajectory  state-vector 
estimation  normally  requires  as  a  minimum,  the  estimation 
of  a  nine  dimensional  state  vector,  three  coordinates  at 
each  of  the  three  levels:  position,  velocity,  and  accelera¬ 
tion.  The  transformations  from  rectangular  to  spherical 
polar  coordinates  (FPS-16  radars)  or  sine-space  coordinates 
(phased-array  radars)  are  non-linear.  However  at  the 
velocity  level  the  transformations  are  achieved  via  the 
Jacobian  matrix  as  the  connection  matrix. 

In  most  statistical  data  processing  procedures  for 
range  instrument  measurements  one  computes  residuals  or 
error  vectors  between  actual  measurements  and  estimated 
measurements  and  actual  states  and  estimated  states. 

These  error  vectors  normally  imply  discrete  differences 
or  approximations  to  differentials  of  vectors.  The  Jacobian 
matrix  provides  the  transformation  between  differentials  of 
rectangular  coordinates  and  spherical -polar  (radar  coordi¬ 
nates)  or  any  other  set  of  generalized  coordinates,  for 
example  sine-space  coordinates  occurring  in  phased-array 
radar  systems.  Conventional  least  squares  polynomial 
filtering,  or  Kalman  filters  in  which  one  linearizes  about  a 
nominal  trajectory  or  linearizes  about  the  current  state 
estimates  can  lead  the  unwary  analyst  into  computational 
problems  of  ill-conditional  matrices,  or  the  scalar  analog 
of  inverting  trigonometric  functions  when  the  variables  are 
near  singularities. 

Quite  often  one  estimates  state  variables  at  the 
rectangular  coordinate  level,  for  example  three  position, 
three  velocity,  and  three  acceleration  coordinates.  At 
some  point  in  the  computations  one  transforms  back  to 


instrument  coordinates  to  compute  instrument  errors,  ’’'he 
computation  of  a  matrix  of  variances  and  covariances  between 
the  variables  implies  two  variance  matrices  one  a  function 
of  the  rectangular  coordinates,  the  second  a  function  of 
the  generalized  coordinates  and  such  that  they  are  related 
through  a  congruent  transformation  on  the  Jacobian  matrix. 

If  one  performs  ;tatistical  analysis  and  computations 
on  data  originating  f;*om  these  instrumentation  systems  he 
should  be  wary  of  ill-conditioned  matrices  or  the  equivalent 
problem  of  solutions  appearing  unstable  or  meaningless. 

The  greatest  insight  into  the  values  of  the  variables 
where  computational  problems  can  occur  comes  from  a  vector- 
matrix  analysis  of  the  geometry. 

It  is  perhaps  well-known  that  singularities  of  FPS-I6 
type  radar  coordinates  occur  near  90°  elevation  angles. 
Perhaps  it  is  not  so  well  known  that  sine-space  (phased- 
array  radars)  singularities  occur  near  45°  angles.  Angle 
measuring  equipment  (AME)  that  measures  direction  cosines 
has  singularities  at  zero  elevation  angles. 

The  mathematics  of  these  relations  are  derived  and 
discussed  below. 
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COORDINATES  AND  JACOBIAN  MATRICES 


The  non-linear  relationships  be  .ween  the  rectangular 
coordinates  (x,  y,  z)  -  or  when  convenient  (xj,  X2>  x3)  - 
and  five  other  sets  of  covering  coordinates  are  developed. 

At  the  position  1'  el  the  functional  relationships  are  as 
always  non-linear.  At  the  differential  level  the  coordinates 
are  connect?'  ;>y  a  Jacobian  matrix,  hence  matrix  analysis 


tools  are  usv.ul, 
theory . 


or  the  tools  of  linear  transformation 


Consider  five  sets  of  covering  coordinates  as  indicated 
by  the  many  variables  shown  in  Figur*:  1. 


FIGURE  1  -  Covering  Coordinates 

The  point  designated  by  the  vector  x  has  three  independent 
coordinates  i=l,2,3,  non-linearly  related  to  three  other 
independent  coordinates  (except  at  singularities)  i=l,2,3 


Xl 

=  X! (q 1 , 

qa» 

qs) 

X2 

=  X2(qi , 

qz. 

q») 

U) 

Xs 

=  X3(qi, 

qz. 

qs) 

for  five  sets  of  q^’s. 
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Designate  the  five  systems  as 


a.  Rectangular  cartesian  coordinates  (x,y,z)  or  (xjf 
x2,  xs). 

b.  Spherical  polar  coordinates  (r,A,  U)  FPS-16  type 
radar.  Optics  type  coordinates  without  range  r. 


c.  Range  and  two  direction  cosines  (r,  cos  ©i,  cos  ©2) 
the  unit  sight  line  vector 


<x 

<xx> 2 


=  <s  --  (slr 


s  J  > 


S3) 


12) 


is  a  function  of  only  two  of  the  three  independent  coordi¬ 
nates  since 


1  *  sf  +  sf  +  sjf  .  (3) 


The  rectangular  coordinates,  Sj_  take  on  many  equivalent  forms, 
for  example 


Si  =  cos  0i  =  CECA  =  etc. 


(4) 


Range  equipment  such  as  AME,  (angle  measuring  electronic 
equipment)  outputs  direction  cosines  (minus  range) . 

d.  Sine-space  coordinates  (r,  sin^i,  sirnjij). 

These  coordinates  occur  in  phased  array  radars.  If  one 
redefines  the  direction  angles  0,  of  Figure  1  in  terms  of 
complements,  that  is 


<J>i  *  90°  -  0* 


(5) 


then 

s  1  =  cos©  1  =  sin<J>i 
sj  =  cos0j  =  sin<J>s 


(6) 

(7) 
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a  number  of  phased-array  radar  studies  appear  to  define 
angles  *  and  8  such  that 


ip  )  =  -a 
<J>5  =  -8 


(8) 


a  number  of  publications  related  to  <*i3  trackers  as  filter 
etc.  exist. 

e.  The  final  set  of  independent  covering  coordinates 
are  range  and  two  of  the  direction  angles  say 


(r ,  i ,  <!),) . 


are 


The  differentials  of  the  scalar  functions  of  Equation  1 


and 


dx1 


\ 

dq. 


(9) 


,2  ,'3x2  ,  \ 

dx  ~  d q/ 


5  % 

If  J<l. 


dx3  = 


\ 


Bx1  = 

.3  q 

(dx1 

\3qi 

dx1 
’  3q  2 

wise 

dx1' 

\ 

.  dq  i 

dx2  1 

=  j 

•  dq  2 

\ 

dx 3 

dq. 

3x 


i\ 


O  A  I 
/ 


(10) 


(11) 
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where  the  Jacobian  matrix  of  (11)  has  the  three  gradient 
vectors  of  (9)  as  its  row-space 

i 9x 1  i 

:  N*  q 


J 

3x3 


/ix± 

q 


ax1; 

"Sqy 


(12) 


■  3x  5 

\3  q 


The  inverse  of  J  exists  when  the  three  gradient  vectors 
are  linearly  independent.  Likewise  when  the  vectors  are 
linearly  dependent  the  matrix  is  singular  and  the  values  of 
the  generalized  coordinates  for  which  singularities  occur 
are  called  poles.  For  example  with  spherical  polar  coordinates 
(r,A,E)  at  the  pole  (thinking  of  the  earth)  elevation  angle  E 
(or  latitude)  equals  90°. 

It  is  well  known  that  the  values  of  the  variables  that 
make  the  determinant  equal  to  zero  are  the  singularities  that 
is 


det  J  =  0.  (13) 

where 

J_1=  adj  J  (det  J)  (14) 


and  adj  J  is  adjoint  of  J  or  transpose  of  matrix  of  co-factors. 
It  is  also  of  interest  to  express  the  inverse  as 

J-1  =  (JTJ)_1  JT  (IS) 


since  the  computational  aspects  of  inversion  is  reduced  to  the 
inversion  of  a  symmetric  matrix  (an  easier  task- in  general). 
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It  is  also  conceptually  worthwhile  to  consider  the  symmetric 
matrix  jTj  and  JJ*  which  in  general  are  not  equal.  These 
"Grammian"  or  "Metric-Matrices"  yield  insight  into  orthogonality, 
obliqueness  etc.  of  row  and  column-space  vectors.  No  further 
structure  properties  are  considered  in  this  report. 

The  notation  Jj^ 

k,  m  =  1,  2,  ...  5 


is  used  to  designated  the  coordinate  system  for  example  J\z 
is  the  Jacobian  matrix  connecting  systems  1  and  2 

/  dx1  \  /  dr 

/  \  \ 

(  dx2  j  *  J12  (  dA  l  .  (16) 

'  dx1  /  \  dE  ' 

In  the  following  tables  some  of  the  Jacobians  are  expressed 
as  functions  of  the  rectangular  coordinates  some  of  the 
generalized  coordinates  and  some  for  both. 

Clearly  if  one  divides  Equation  (16)  by  At  and  passes  to 
the  limit  by  velocity  relations  result,  that  is 


'  k1  N 

qi  ' 

J  \  f 
s 

\ 

•  =  J  ( 

q  2  ; 

(17) 

\  \ 

/ 

« 

\  x’ 

^  qs 

where  the  example  of  (16)  yields 

q  =  (f.  A,  £).  (18) 
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Variance  matrix  relations  are  also  obtained  by  trans¬ 
posing  Equation  (11) 

dx  a  <^dqJx  (39) 

and  computing  the  dyadic  product 

d£><$lx  55  Jd<J>sdqJ^  (20) 

In  a-priori  flight  testing  one  has  a  sequence  j  of  past 
tests  and 

,  Jiaax 

£(dx><dx)  =  4 -  I  dx><dx  =  l  (21) 

Jmax  j=l  3x3 

or 

Lex  -  JlqqJT  (22) 

where  e  is  expectation  operator. 

Sample  variance  estimates  instead  of  past  or  population 
sequences  yield  similar  equation  forms.  The  Jacobians  are 
evaluated  about  a  nominal,  or  known  estimated  point.  Equation 

(22)  expresses  the  well  known  variance  connections  via  con¬ 
gruent  transformations. 

In  order  to  further  simplify  the  typing  and  the  results 
observe  that  the  unit  magnitude  "sight- line"  vector  is 


‘ (Si’ S2’ 


(23) 


where 


<s§>  =  1 


(24) 


or 


and 


s2 


Ss 


•  -rr- 


si 


_  2 
i>  S 


nJ  l-sf 


_  * 
s  2 


(25) 


(26) 


The  rectangular  coordinates  called  "direction  cosines'*  of 
the  sight-line  vector  has  many  equivalent  forms,  some  of 
which  are 


/S"v  ^ CECA \  ^c0i  ^  /  cGk 
|  s2  N  CEsA  j=|  sa  -  c02  | 

^sE  S  VC0S/  Vss  J 


(27) 


where  the  square  root  terms  of  Equation  (25)  and  (26)  aie  used. 
The  following  matrices  are  developed 


d  1  2 
J  1  S 
J 1  l) 
Jis 


Ja  s 
J  2  «• 

J2s 


J 


3  k 


1  «i  5 


3  5 


Consider  the  singularities  (JiO  in  transforming  from 
sine-space  to  rectangular  coordinates.  The  singularities 
occur  when 


S2$s  *  i 


or  vhei 
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v>  rt 


4>i  +  4>3  =  99c 


hus  we  see  a  continum  of  singularities  as  opposed  to  a 
ingle  point  (£=90“)  for  conventional  radar  or  spherical 
polar  coordinates. 

One  trajectory  region  to  watch  out  for  ill-conditioned 
matrices  is  when  <j)i=<f>2*45° . 

The  analyst  can  study  the  singularity  regions  as 
applicable  to  his  particular  problem.  The  position  vectors 
between  station  sites  are  not  included  in  this  report  but 
will  be  developed  in  a  later  report. 
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CONCLUSIONS 


The  geometi ieally  related  computational  problems  pre¬ 
sented  here  are  only  part  of  the  relations  occurring  in 
range  data  processing.  To  be  sure,  the  position  vectors 
between  sites  and  the  relations  between  generalised 
coordinates  measured  at  the  two  different  sites  must  be 
considered.  Error  propagation  studies  for  a  total 
instrumentation  system  will  be  developed  and  draw-upon 
seme  of  the  relations  herein. 
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